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We investigate properties of r-mode instability in slowly rotating relativistic polytropes. Inside the 
star slow rotation and low frequency formalism that was mainly developed by Kojima is employed to 
study axial oscillations restored by Coriolis force. At the stellar surface, in order to take account of 
gravitational radiation reaction effect, we use a near-zone boundary condition instead of the usually 
imposed boundary condition for asymptotically flat spacetime. Due to the boundary condition, 
complex frequencies whose imaginary part represents secular instability are obtained for discrete 
r-mode oscillations in some polytropic models. It is found that such discrete r-mode solutions can 
be obtained only for some restricted polytropic models. Basic properties of the solutions are similar 
to those obtained by imposing the boundary condition for asymptotically flat spacetime. Our results 
suggest that existence of a continuous part of spectrum cannot be avoided even when its frequency 
becomes complex due to the emission of gravitational radiation. 



I. INTRODUCTION 



Andersson jl[ and Friedman and Morsink Q discovered that all r-modes, which are quasi-toroidal modes mainly 
restored by Coriolis force, in all rotating stars become unstable due to the gravitational radiation reaction if other 
dissipative processes are not considered. This instability is clearly understood by the so called CFS mechanism |5]-[|. 
As shown by Lindblom, Owen, and Morsink || first, this instability still strongly affects on stability of typical neutron 
star models even if viscous dissipation of neutron star matter, which tends to stabilize the CFS instability, is taken 
into account. Since then a lot of studies on oscillation modes restored by Coriolis force in rotating stars have been 
done to prove their possible importance in astrophysics (for recent review, see, e.g., Rcfs 00). 

Influence of the r-mode on stability of rotating neutron stars is one of the most important and interesting phe- 
nomenon in astrophysics. In oscillations of neutron stars, relativistic effect must be important because such stars 
are sufficiently compact. But most studies have been done within the framework of Newtonian gravity so far, al- 
though our understandings of r-modes have been improved by those investigations. As for r-modes studied within 
the framework of general relativity, Kojima jll] derived master equations for r-mode oscillation in the lowest order 
slow rotation approximation, and then he found possible existence of a continuous part of spectrum in his equations. 
Beyer and Kokkotas generally verified the existence of a continuous part of spectrum in Kojima's equation. Ko- 
jima's formalism was developed to include higher order rotational effects by Kojima and Hosonuma [ ^3|Jl4| . Lockitch, 
Andersson, and Friedman obtained the discrete r-mode solutions in uniform density stars as well as a continuous 
part of spectrum by solving Kojima's equation. Recently Yoshida and Ruoff and Kokkotas Jl?]] discussed that 
such discrete r-mode solutions are not simply allowed to appear in compressible stellar models. Their results showed 
that for typical neutron star models, Kojima's equations do not have such a discrete r-mode solution. 

These recent developments of understanding of relativistic r-modes have shown that basic properties of r-mode 
oscillations in relativistic stars are significantly different from those in Newtonian stars. As for non-isentropic stars, 
most previous studies have shown the existence of a continuous part of spectrum. This is a great contrast with 
Newtonian case. For Newtonian cases, there are discrete mode solutions and no continuous parts of spectrum for 
r-modes in all uniformly rotating stars as long as their rotation velocity is small enough. It is not however likely 
to occur such drastic change in the behavior of the solutions due to the inclusion of even a little relativistic effect. 
Therefore, most authors have considered that a continuous part of spectrum does not appear if some effects that were 
omitted in previous studies are taken into account. One of such effects is a dissipation effect due to gravitational 
radiation reaction. In most studies on r-modes, a slow motion approximation has been employed in their analysis 
because of slow rotation approximation. A slow motion approximation changes a wave type equation into a Laplace 
type equation. Thus, Kojima's equations do not have solutions with wave character, and then the frequency is real 
number if asymptotically flat spacetime is assumed. In other words, influence of gravitational radiation on relativistic 
r-modes has not been taken into account so far. In this case, Kojima's equation becomes that of singular eigenvalue 
problem for some frequency range, and hence has a continuous part of spectrum. As suggested by Lockitch et al. 
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|L5| (see also, Rcf. |l2| ), however, Kojima's equation may become that of regular eigenvalue problem if frequency has 
non-zero imaginary part. Therefore, It is hoped that the existence of a continuous part of spectrum might be avoided 
if the effect of the emission of the gravitational radiation is considered. 

In this paper, accordingly, we will attempt to include the lowest order effect of gravitational radiation reaction into 
Kojima's formalism. Because Kojima's equation is not wave type equation, as mentioned before, we can not obtain 
information about gravitational wave from the equation at all. In order to include the gravitational radiation reaction 
effect into Kojima's equation, we will employ a near-zone boundary condition instead of usually used boundary 
condition for asymptotically flat spacetime. This boundary condition was introduced by Thorne | fl8[ to include the 
effect of gravitational radiation reaction on polar pulsations in a Newtonian star. We start, in §2, with the description 
of our method of solution. A near-zone boundary condition is introduced there in order to take account of the 
gravitational radiation reaction. In §3, we show the properties of the r-mode solutions obtained by imposing the 
boundary condition. §4 is devoted for discussions and conclusions. Throughout this paper we will use units in which 
c = G = 1, where c and G denote the velocity of light and the gravitational constant, respectively. 



II. METHOD OF SOLUTIONS 



We consider slowly rotating relativistic stars with a uniform angular velocity 17, where we take account of the first 
order rotational effect in f2. The geometry around the equilibrium stars can be described by the following line clement 
(see, e.g. Ref. ||): 

ds 2 = -e 2l,{r) dt 2 + e 2Mr) dr 2 + r 2 d0 2 + r 2 sin 2 Odip 2 - 2oj(r) r 2 sin 2 6dtd<p . (2.1) 

Throughout this paper, the polytropic equation of state is assumed: 

p = Kp 1+ ^, (2.2) 

where p and p denote the pressure, mass-energy density, respectively. Here N and K are constants. 

We consider oscillation modes in rotating relativistic stars such that the eigenfunctions are stationary and are 
composed of only one axial parity component in the limit Q — > 0. Subclass of those modes should be a relativistic 
counterpart of r-modes, which are able to oscillate in all slowly rotating Newtonian fluid stars. According to Lockitch 
ct al. |l5|] , such modes are allowed to exist only if the star has non-isentropic structures. Therefore, we assume stars to 
be non-isentropic, although the effects due to deviation from isentropic structure on oscillation modes do not appear 
in the first order in fi. According to the formalism by Lockitch et al. Jl5|] (see, also, Ref. [0), let us write down the 
pulsation equation for relativistic r-modes with accuracy up to the first order in J7. The metric perturbation 8g a p 
and the Eulerian changes of the fluid velocity Su a that do not vanish in the limit O — > are given as 

d v Y lm (6,ip) \ - t 



(5g te ,Sg tv ) = ih ,i(r) - y v „ ,8m6d e Y lm (6,<p) e lat , (2.3) 



r z \ smfl smtf / 

where Yi m (6, <p) are the usual spherical harmonic functions, and a denotes oscillation frequency measured in the 
inertial frame at the spatial infinity. All other perturbed quantities become higher order in £1. Note that this form 
of eigenfunctions is the same as that for zero- frequency modes in a spherical non-isentropic star, because r-modes 
become zero-frequency ones in the limit O — >• |20[| . The metric perturbation /lo.z obeys a second order ordinary 
differential equation, 



Di m (r;a) 
where 



-a d ( _ v _ x dha,i\ (Kl + 1) , -2 + 2e- 2A 



€ dr ( e * ) - I + + 8lT{P + P) ) k0 > 1 



+ 16w(p + p)h ,i=0, (2.5) 



D lm (r;a) = l- , 2.6 

1(1 + 1)<7 



Here, we have introduced the effective rotation angular velocity of fluid, 
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ui = D,-oj, (2.7) 

and the corotating oscillation frequency, 

a = a + mtt. (2.8) 
The velocity perturbation of fluids Ui is determined from the function hoj through the following algebraic relation, 

2mu) 



1 



l{l+\)a 



Ui + h ,i = . (2.9) 



Equ ations (2.5) and ( [2.9] ) are our basic equations, which were derived by Kojima |ll|] first. Note that the equations 



(2.5) and (2.9) lose th eir meaning in the I = case, because there are no axial modes with I = 0. 



Because equations (2.5) are second order ordinary differential equations, two boundary conditions are required to 
determine solutions uniquely. From regularity of physical quantities at r = the function ho t i must vanish at the 
center of a star. This condition is explicitly given as 

r do£_, l + 1 s hQi= Q asr ^ (2.io) 
dr 

Outside the star equations ( |2.5| ) have general solutions as follows: 

oo oo 

h ,i(r) =AY,a k r- l - k + S^6 fc r i+1+fc , (2.11) 

fc=0 fc=0 

where A and B are arbitrary constants. Here au and bk are constants determined from recurrence relations although 
explicit expressions are omitted here. Note that and do not depend on frequency a. In most previous studies the 
condition B — have been chosen as a boundary condition at the spatial infinity because spacetime must be regular 
everywhere. In this paper we call this condition "proper boundary condition" . This condition means the spacetime is 
asymptotically flat. Therefore, solutions satisfying the condition B — can not describe any gravitational radiation 
emission from a star at all. In order to include effect of gravitational radiation emission in the solutions, we must not 
require to demand the zero for B. This was first shown and used in the study for post-Newtonian approximation by 
Thorne @. 

Next, let us consider the boundary conditions for quasi-normal mode solutions, according to similar consideration 
to that by Thorne Jl^] (see, also, Ref. pi]). In the derivation of our basic equations slow motion approximation is 
necessarily used as well as slow rotation one because the frequency of oscillation restored by Coriolis force is the same 
order of the stellar rotation frequency. This slow motion approximation must nicely work near the star as long as low 
frequency oscillations are considered. Therefore, higher order time derivative of perturbed quantities are neglected in 
governing equations ( |2.5| ). If an oscillating star emits gravitational radiation, however, some of such omitted terms 
must become important in the radiation zone because a term such as err becomes dominant among all terms in the 
governing equations. Besides, rotational effects due to the stellar rotation fall off faster than 1/r 2 as r — > oo. Thus, 
such slow motion approximation is not good far from the star, if gravitational wave are radiated from the stars. 
Regge- Wheeler equations with correction terms due to the stellar rotation, in fact, govern the axial perturbations 
sufficiently far from the star even when low frequency oscillations like r-modes are induced. Since err > 1 in radiation 
zone, Regge- Wheeler equations can be approximately written as 



r 



2 «(r) : ,„ , , „(M 



dr 2 



+ a 2 r z Xi(r) - 1(1 + l)Xi(r) + O — = , (2.12) 



where Xi(r) are Regge- Wheeler's functions. Here M denotes the gravitational mass of the star. The metric functions 
/to,; are determined from the functions Xi by the equation: 



dr \ r J 

(see, e.g. Refs p2|,p3[ ). Here we have considered the limiting case when M/r <C 1 for simplicity. The general solutions 
to equations ( |2.12] )can be given analytically: 

Xi(or) =ar(Cji(ar)+Dm{ar)), (2.14) 
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where ji and n; are spherical Bessel functions and C and D are arbitrary constants. The asymptotic forms in the 
radiation zone, that is, when err > 1, are given as 



Xi(ar) ~ C cos 



1)tt 



I? sin 



or — — (I + 1)tt 



(2.15) 



Now we are interested in quasi-normal modes of a star. Thus, the no incoming radiation conditions are required for 
metric perturbations. From equation (2.15) it is found that this condition becomes the relation D — —iC. For this 
choice of constants the asymptotic solutions reduce to the form as: 



X x (or)e iat ~ C exp 



io(t -r) + -(l+ 1)tt 



as or 



(2.16) 



Since the frequencies of the oscillations restored by the Coriolis force are proportional to the rotational frequency Q, 
o 2 R 2 ~ Q 2 R 2 = t 2 M/R, where e denotes a small parameter for stellar rotation defined as e = Q/ (Mj i? 3 ) 1 / 2 . Here, 
R is the radius of the star. Thus, the surface of the star is approximately in the near zone, that is, or <C 1, if the 
stellar rotation is sufficiently slow or the stellar gravity is sufficiently weak. In this approximation, solutions (2.14) 
can be written as 



Xi(or) 



■iC 



(gj - 1) 
(or) 1 



l + i 



(or) 



21+1 



(2Z-1)!!(2Z + 1)!! 



(2.17) 



where the constraint D = —iC for the no incoming radiation has been used. The approximate solutions above may 
be valid near the stellar surface because M/r < 1 is well satisfied at the stellar surface and outside the star for typical 
neutron star models. In the near zone, thus, the expressions for metric perturbations ft-o,; outside the star can be 



given by equations ( 2.13 ) and ( 2.17 ) as follows: 



ho,. 



l + i 



(l + 2)(or) 



21 + 1 



(i-l)(2i + l)[(2i-l) 



1112 



(2.18) 



where C is an arbitrary constant. As the outer boundary condition we require this solution to connect smoothly to 
interior solution obtained by solving equation ( |2.5| ) at the stellar surface. Thus, the boundary condition is explicitly 
given as 



1 



(l + 2)(oR) 



21+1 



(l-l)(2l + l)[(2l-l)»}\ 
. (l + l)(l + 2)(oR) 2l+1 
" Z (Z - 1)(2Z + 1)[(2Z - 1)!!] 



r^i(R- x ) 
dr 

h ,i(R-x) = 0. 



0. 



(2.19) 



where we need actual values of rotation frequency to obtain eigensolutions. In this paper we assume the value of 
as fl = (71-p) 1 / 2 , where p is the average density defined by p = M/(47ri? 3 /3). This f2 is an approximate value for 
maximum rotation frequency to be po ssible to settle down uniformly rotating stars in hydrostatic equilibrium. In 
this paper, we call the condition (2.19) "near zone boundary condition". This method is a crude version of matched 
asymptotic expansions. If we consider a non-rotating limit, that is, o = limit, the boundary condition (2.19) becomes 
approximation of proper boundary conditions in which only the lowest order terms in M/R are included. As shown 
by Lindblom et al. pl| , the boundary condition similar to that by using asymptotic solution (2.18) can give a good 
approximate value of eigenfrequency even for /-mode oscillation although /-mode is not a low frequency oscillation 
mode. 



III. NUMERICAL RESULTS 



As shown in the previous studies |1 l| , |l2| , [l5f , we should distinguish two cases in treating equation (2J3) when the 
proper boundary condition is imp osed. One case is regular eigenvalue problem and the other singular one. For regular 
eigenvalue problem equation (2.5) may have discrete mode frequency in a range, 



2muj(R) _ 2mw(oo) 2mf2 

W+-T) <a - 



1(1 + 1) i(i + i) 



(3.1) 



4 



On the other hand, equation (2.5) becomes singular eigenvalue problem if a is in a region, 



2m^(0) _ 2mg(g) 

7(T+ij- - i+iy ( } 

Notice that the range ( |3.2[) is a continuous part of spectrum of equation (|2.5[) . As pointed out by Lockitch et al. 



(see also, Ref. |12|) equation (2.5) becomes that of regular eigenvalue problem when the corresponding frequency has 
non-zero imaginary part. Thus frequency ranges above may not have clear mathematical meanings when frequency 
has non-zero imaginary part. According to frequency ranges above, however, we will divide eigensolutions into three 
classes: The first and the second class solutions are characterized by their real part of frequency satisfying inequality 



( |3.ip and (3.2), respectively. The third class composes of a compensative set of the first and the second class. 

First of all, we concentrate our attention to r-mode solutions with frequency whose real part is in a range ( ^l| ) . We 
compute frequencies of mode solutions for several polytropic stellar models. In the present study, only the fundamental 
r-modes, whose eigenfunction U m has no node in radial direction except at stellar center, are obtained. This is similar 
to that in studies for the proper boundary condition case [ ]l6| , p"7| |. In Figures 1 and 2, real and imaginary parts of 
scaled eigenfrequencies n = a/fl of r-modes are, respectively, given as functions of M/R. Eigenfrequencies for stars 
with four different polytropic indices, N — 0, 0.5, 0.75, and 1, are shown in panels in both figures, respectively. 
Only frequency curves for the modes with I = m = 2 are depicted along a relativistic factor M/R because they are 
considered to be the most important modes for r-mode instability. 

Real parts of frequency illustrated in Figure 1 are in good agreement with Figure 1 of Ref. jl6| , in which frequency 
are obtained by imposing the proper boundary condition for asymptotically flat spacetime. The relative differences 
are less than 0.1%. This shows that our approximation works nicely and higher order effect of M/r on the outer 
boundary condition is not so important for the determination of real parts of frequency. We also find that frequency 
curves in Figure 1 are terminated at some value of M/R beyond which equilibrium states can still exist. Here, the 
maximum values of M/R for polytropic equilibrium stars having N = 0, 0.5, 0.75, and 1.0 are given by 4/9, 0.385, 
0.349, and 0.312, respectively. It is also found that length of frequency curves tend to be short as a polytrope index 
N increases. This feature is similar to that for the case where the proper boundary condition is used. And those 
terminal points of frequency curves appear at almost the same values of relativistic factors as those for the proper 
boundary condition case (see Ref. Jl(|). Beyond the value of the relativistic factors corresponding to those terminal 
points, we can obtain a lot of eigensolutions with a singular eigenfunction but not with a regular one. Furthermore, 



the real part of the corresponding frequency belongs to a range (3.2) but not a range (3.1). 

From Figure 2, we can see that r-modes obtained in this study are all unstable. It is also found that curves for 
imaginary parts of k have one relative minimum near the terminal point of frequency curves. Those minimum values 
are given as lm(n(M/R = 0.425)) = -2.9 x 10~ 3 , lm(n(M/R = 0.344)) = -6.4 x 10~ 4 , lm(n(M/R = 0.258)) = 
-1.3 x 10~ 4 , and \xq.{k{M/R = 0.095)) = -2.2 x 10~ 6 , for stars having N = 0, 0.5, 0.75, and 1.0, respectively 
Values of Im(«;) also approach zero as relativistic factor M/R is getting closer to a value corresponding to that 
for the terminal point of frequency curves. Those behaviors of Im(K) can be understood from the distribution of 
eigenfunctions Ui(r) because Im(K) is approximately proportional to the square of the current multipole moment. In 
Figure 3 and 4 distributions of the eigenfunction U m are shown for N — 0.5 polytropic models having M/R = 0.1 
and M/R = 0.37, respectively. As seen from these figures the motion of perturbed fluid elements is strongly confined 
near the stellar surface when the mode frequency is getting closer to that of the terminal frequency, which satisfies 
(7 w 2muj(R) / (1(1 + 1)). We can easily understand this behavior from equation ([2.9]). Consequently the values of 
current multipole moments of such modes may become small when the value of relativistic factor increases. On 
the other hand, the efficiency of the gravitational radiation emission becomes good as increase of the value of the 
relativistic factor. Due to both effects above a relative minimum of imaginary part of frequency may appear. We 
should notice that in N = case, a value of Iul(k) do not approach zero even when M/R ~ 4/9, as we can see from 
Figure 2. The reason is that the eigenfunctions U m does not have so strong peak at the stellar surface even when 
M/R ~ 4/9 because a relation a w 2m0j(R) / (1(1 + 1)) is not satisfied in this case. 

Next let us estimate the instability timescale for the gravitational radiation driven instability of r-mode solutions 
with frequency in a range ([O]). Here a typical neutron star model whose mass and radius are respectively 1.4M Q and 
12.57 km is considered for N = 0, 0.5, and 0.75 polytropic models. When the star rotates with the angular frequency 
n = (3/4 x GM/R?) 1 / 2 = 8377 s"\ growing timescales Tj of the r-mode instability are given as Tj = 1.29 s, 2.04 s, 
and 2.97 s for N = 0, 0.5, and 0.75 models, respectively. Note that for N = 1 case the growing timescale cannot 
be estimated because we can not find discrete r-mode solution for that model. These timescale are similar to those 
obtained from Newtonian estimate, (see, e.g. Ref. 0) 



When we concentrate our attention to solutions whose frequency has real part in a frequency range (3.2), we obtain 
a large number of solutions whose eigenfunction has singular behavior in its real part at a radius determined by a 
solution of Di m (r; Re(a)) = 0. Besides, those solutions have severe truncation error due to a finite difference method. 
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Similar behavior of solutions appears in oscillations of differentially rotating disks, (see, for example, Ref. (24]) As 
discussed by Schutz and Verdaguer |24| , this is considered a sign of the existence of continuous parts of spectrum, 
although the existence has to be proved by other mathematical techniques because exact continuous spectrum is 
never obtained from a simple numerical analysis. In the present case, we are sure that the appearance of a continuous 
spectrum is plausible because the imaginary parts of frequency are too small to change drastically the character 
of the solutions derived from the proper boundary condition. Thus, our numerical results suggest the existence of 
a continuous part of spectrum in Kojima's equa tion even when their frequency becomes complex number. As for 
regular solutions with frequency in a range (3.2), we cannot obtain such a solution at all. Finally we consider the 
third class of solutions, whose real part of frequency is neither in a range (3.1) nor (3.2). In this region of real part 
of frequency, we cannot obtain any solutions at all. 



IV. DISCUSSION AND CONCLUSION 



In this paper, we have investigated the properties of r-mode instability in slowly rotating relativistic polytropes. 
Inside the star slow rotation and the low frequency formalism that was mainly developed by Kojima [Tl| and Lockitch 
et al. is employed to study axial oscillations restored by Coriolis force. At the stellar surface, in order to take 
account of gravitational radiation reaction effect, we use a near-zone boundary condition, which was devised by 
Thorne jl8| and recently developed for relativistic pulsations by Lindblom et al. ]2l]] , instead of the usually imposed 
boundary condition for asymptotically flat spacetime. Due to the boundary condition, complex frequencies whose 
imaginary part represents secular instability are obtained for r-mode oscillations in some polytropic models. It is found 
that such discrete r-mode sol utio ns can be obtained only for some restricted polytropic models. Basic properties of 
mode solutions of equation (2J3) that obtained in this study is similar to those with the boundary condition for 
asymptotically flat spacetime although its frequency becomes complex because of the near-zone boundary condition. 



As suggested by Lockitch et al. |15| (see, also Refs 12,lfl,h7]]), when an eigenfrequency becomes complex number, 



which expresses the damping of the oscillation due to the energy dissipation such as the gravitational radiation 
emis sion, there is a possibility that the existence of the continuous part of spectrum in the eigenfrequency of equation 
(2.5) is avoided. In this study, we consider the complex frequency corresponding to the quasi-normal mode as mode 
solutions of equation (2.5) by imposing the lowest order near-zone boundary condition. Our numerical results suggest 
the existence of a continuous part of spectrum in Kojima's equation even when the frequency is allowed to be a 
complex number. However, we still think that the existence of a continuous part of spectrum in axial oscillations 
restored by Coriolis force is not plausible because such property does not appear in Newtonian r-modes. The existence 
of a continuous part of spectrum in this study might be artifact due to the approximation because our treatment is 
the lowest order approximation. In other words, inclusion of full effect of gravitational radiation emission might avoid 
the existence of a continuous part of spectrum. Another possibility to prevent the appearance of continuous spectrum 
might be to avoid a singular reduction in the order of the equation. As Kojima and Hosonuma |fl4| showed, basic 
equations for r-mode oscillations become a fourth order ordinary differential equation for metric perturbation ho; 
when rotational effect up to the third order of ^/(M/i? 3 ) 1 / 2 is consistently considered. In this equation, extra two 
degrees of freedom of solutions may be used to avoid singular behavior of eigenfunction. Due to these extra boundary 
conditions, all eigenfrequency may become discrete. Verification of those possibilities remains future studies. 
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FIG. 1. Real parts of scaled frequencies k = a/fl of the r-modes with I = m — 2 plotted as functions of M/R. In each panel, 
the frequencies of modes for polytropic models with N = 0, 0.5, 0.75, and 1 are respectively shown. The labels indicating 
their polytropic indices N are attached in corresponding panels. Vertical dotted lines show the maximum values of M/R for 
equilibrium states: M/R = 0.444 for N = 0, M/R = 0.385 for N = 0.5, M/R = 0.349 for N = 0.75, and M/R = 0.312 for 
N = 1.0. 
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FIG. 2. The same as Figure 1 but for imaginary parts of frequency. 
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